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Abstract 

In this paper, we consider a nonlinear inverse problem occuring in nu¬ 
clear science. Gamma rays randomly hit a semiconductor detector which 
produces an impulse response of electric current. Because the sampling 
period of the measured current is larger than the mean interarrival time 
of photons, the impulse responses associated to different gamma rays can 
overlap: this phenomenon is known as pileup. In this work, it is as¬ 
sumed that the impulse response is an exponentially decaying function. 
We propose a novel method to infer the distribution of gamma photon en¬ 
ergies from the indirect measurements obtained from the detector. This 
technique is based on a formula linking the characteristic function of the 
photon density to a function involving the characteristic function and its 
derivative of the observations. We establish that our estimator converges 
to the mark density in uniform norm at a polynomial rate. A limited 
Monte-Carlo experiment is provided to support our findings. 


1 Introduction 

In this paper, we consider a nonlinear inverse problem arising in nuclear science: 
neutron transport or gamma spectroscopy. For the latter, a radioactive source, 
for instance an excited nucleus, randomly emits gamma photons according to 
a homogeneous Poisson point process. These high frequency radiations can be 
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associated to high energy photons which interact with matter via three phe¬ 
nomena : the photoelectric absorption, the Compton scattering and the pair 
production (further details can be found in [14]). When photons interact with 
the semiconductor detector (usually High-Purity Germanium (HPGe) detectors) 
arranged between two electrodes, a number of electron-holes pairs proportional 
to the photon transferred energy is created. Accordingly, the electrodes gen¬ 
erate an electric current called impulse response whenever the detector is hit 
by a particle, with an amplitude corresponding to the transferred energy. In 
this context, a feature of interest is the distribution of this energy. Indeed, it 
can be compared to known spectra in order to identify the composition of the 
nuclear source. In practice, the electric current is not continuously observed but 
the sampling rate is typically smaller than the mean inter-arrival time of two 
photons. Therefore, there is a high probability that several photons are emitted 
between two measurements so that the energy deposited is superimposed in the 
detector, a phenomenon called pile-up. Because of the pile-up, it is impossible to 
establish a one-to-one correspondence between a gamma ray and the associated 
deposited energy. 

This inverse problem can be modeled as follows. The electric current gen¬ 
erated in the detector is given by a stationary shot-noise process X = 
defined by: 

Xt= ^ Ykh{t-Tk), (1) 

k:Tk 

where h is the (causal) impulse response of the detector and 

(SN-1) Yf, is a Poisson point process with times Tk G M. arriving homoge¬ 

neously with intensity A > 0 and independent i.i.d. marks Yfc € R having 
a probability density function (p.d.f.) 6 and cumulative distribution function 
(c.d.f) F. 

We wish to estimate the density 6 from a regular observation sample Xi ,..., 
of the shot noise (1). Note that the sampling rate is set to 1 without mean¬ 
ingful loss of generality. If a different sampling rate is used, e.g. we observe 
Xs, ..., XnS for some <5 1, it amounts to change A and to scale h accordingly. 

The process (1) is well defined whenever the following condition holds on 
the impulse response h and the density 6 

J min(l, \y h{s)\)9{y)dyds < oo . (2) 


As shown in [12], this condition is also necessary. Moreover, the marginal 
distribution of X belongs to the class of infinitely divisible (ID) distributions 
and has Levy measure v satisfying, for all Borel sets B in M \ {0}, 


A /’°° 

v{B)=X f‘{h{s)YiGB) 

Jo 


ds 


(3) 


The ID property of the marginal distribution shows that this estimation 
problem is closely related to the estimation of the Levy measure v. This property 
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strongly suggests to use estimators of the Levy triplet, see for instance [16] and 
[8]. However, up to our best knowledge, these estimators use the increments of 
the corresponding Levy process which are i.i.d. and they assume a finite Levy 
Khintchine measure. In contrast, the observations are not independent and the 
Levy measure of the process is infinite since from (3), we have that 

poo 

i/(M) = A / P {h{s)Yi £ R) ds = oo. (4) 

Jo 

In order to tackle this estimation problem, we then propose to bypass the es¬ 
timation of ly and directly retrieve the density 6 of the marks distribution F 
from the empirical characteristic function of the measurements. Coarsely speak¬ 
ing, using (3) , the Levy-Khintchine representation provides an expression of the 
characteristic function (px of the marginal distribution as a functional of 9. The 
estimator is built upon replacing px by its empirical version and inverting the 
mapping 9 i—>■ px ■ A more standard marginal-based approach would be to rely 
on the p.d.f. of X. However, the density of Xq is intractable, which precludes the 
use of a likelihood inference method. Consequently, although shot-noise models 
are widespread in applications (for example, such models were used to model 
internet traffic [1], river streamflows [4], spikes in neuroscience ([11],[10]) and in 
signal processing ([19] , [20])), theoretical results on the statistical inference of 
shot-noise appear to be limited. Recently, Xiao and al. ([24]) provide consistent 
and asymptotically normal estimators for parametric shot-noise processes with 
specific impulse responses. 

In this contribution, we consider the particular case given by the following 
assumption. 

(SN-2) The impulse response h is an exponential function with decreasing rate a > 0 : 

h{t) = e““* iR^(t) . 

Under (SN-2), the process (Xt)tgR is usually called an exponential shot-noise. 
In this case. Condition (2) becomes 

E[log+(|yi|)] <oo. (5) 

Under (SN-2), the process (Xt)t>o can alternatively be introduced by consider¬ 
ing the following stochastic differential equation (SDE) : 

dXt = —aXtdt d- dLt , Xg = x G R (6) 

where L = {Lt)t>o is a Levy process defined as the compound Poisson process 

Nxit) 

Lt=Y,Yk with XA(t) = ^lT.<t, (7) 

k—1 k 

where (Tfc,Tfc)fc>o satisfies (SN-1). The solution to the equation (6) is called a 
Ornstein-Uhlenbeck(O-U) process ([18, Chapter 17]) driven by L with initial 
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condition Xq = x and rate a. Note that L defined by (7) has Levy measure XF. 
Thus, by [18, Theorem 17.5], this Markov process admits a unique stationary 
version if (5) holds, and this stationary solution corresponds to the shot-noise 
process (1). 

In recent works, [3, Brockwell, Schlemm] exploit the integrated version of 
(6) to recover the Levy process L and show that the increments of L can be 
represented as: 


pnh 

Lnh — L(^n-l)h = XCnh ~ X(^n-l)h + Q: / Xjds. 

J {n—l)h 

These quantities are only well estimated for high frequency observations so that 
we cannot rely on this method in our regular sampling scheme. 

To the best of our knowledge, the paper that best hts our setting is [13]. The 
authors propose a nonparametric estimation procedure from a low frequency 
sample of a stationary 0-U process which exploits the self decomposability prop¬ 
erty of the marginal distribution. The authors construct an estimator of the so 
called canonical function k defined by: 

i'{dx) = -dx. 

X 

The two main additional assumptions are that k is decreasing on (0,oo) and v 
satisfies the integrability condition /g°°(l A x)v{dx) < oo. In our setting (i.e. 
when specifying the Levy process to be the compound Poisson process defined 
in (7)), it is easily shown that these conditions hold and the canonical function 
and the cumulative distribution of the marks are related by the equation: 

k{x) = AP(Fo >x)=X{l- F(x)). 

In this article, we introduce an estimator of 0 based on the empirical char¬ 
acteristic function and a Fourier inversion formula. This algorithm is numeri¬ 
cally efficient, being able to handle large datasets typically used in high-energy 
physics. Secondly, we establish an upper bound of the rate of convergence of our 
estimator which is uniform over a smoothness class of functions for the density 

e. 

The paper is organized as follows. In Section 2, we introduce some prelim¬ 
inaries on the characteristic function of an exponential shot-noise process and 
provide both the inversion formula and the estimator of the density 9. In par¬ 
ticular, we derive an upper bound of convergence for our estimator over a broad 
class of densities under the assumption that A/a is known. In Section 3, we 
present in details the algorithm used to perform the density estimation and il¬ 
lustrate our findings with a limited Monte-Carlo experiment. Section 4 provides 
error bounds for the empirical characteristic function based on discrete-time ob¬ 
servations and exploit the /3-mixing structure of the process. Finally, Section 5 
is devoted to the proofs of the various theorems. 
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2 Main result 


2.1 Inversion formula 

As mentioned in the introduction, it is difficult to derive the probability density 
function of the stationary shot-noise unless the marks are distributed according 
to an exponential random variable and the impulse response is an exponential 
function. In this case, it turns out that the marginal distribution of the shot- 
noise is Gamma-distributed (the reader can refer to [2] for details). In all other 
cases, we can only compute the characteristic function of the marginal distri¬ 
bution of the stationary version of the shot-noise when treating it as a filtered 
point process (see for example [17] for details). We have for every real u: 

ip[u) = E[e™^«] = exp j - l)duF(d 2 /)^ . (8) 

From (8), the characteristic function of Xq can be expressed as follows: 

‘fXo {u) = exp XKh{uy)F{dy)^ . (9) 

where Kh, the kernel associated to h is given by: 

r-l-oo 

Kh{x) = / (e“'‘(") - l)d^;. 

Jo 

Note that if h is integrable, then is well defined since, for any real x, 
— Ijds < |a:|/ q°° |h(s)|ds. Moreover, if h is integrable, then Kh is 
a C^(R, C) function whose derivative is bounded and equal to: 



Furthermore, if E [lyol] < oo, then the characteristic function of Xq is differen¬ 
tiable and we have: 


f'xoiu) = X y:>Xo{u) [ yK'h{uy)F{dy). (10) 

Jr 

Under (SN-2), the kernel Kh takes the form 

\ /■“ e‘® - 1 

Kh{u)= (e*“^ -l)dv= -ds. (11) 

Jo ^ J Jo “S 

With (10), we obtain that 

^'xoiu) = ^Xo{u)—i^Yoiu) - 1) ■ ( 12 ) 

^ n/oi 
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Since the marginal distribution of X is infinitely divisible, we have by [18, 
Lemma 7.5.] that ipx(u) does not vanish. If in addition ipy is integrable, (12) 
provides a way to recover 9 knowing a/X, namely, for all a; G K, 


^(^) = ^ / {u)du [ e-‘"“ (l + ^ du . (13) 

27r Jr Stt Jr \ A ipxo\u)J 


This relation shows that the estimation problem of the p.d.f. 9 is directly related 
to the estimation of the second characteristic function. 


Remark 2.1. We assume in the following that the ratio a/A appearing in 
the inversion formula (13) is a known constant, as it typically depends on the 
measurement device. Interestingly, however, an estimator of this constant can 
be derived from [12, Theorem 1], where it is shown that the marginal distribution 
G of the stationary shot-noise is regularly varying at 0 with index A/a, i.e. : 

G{x) ~ x^^°‘L{x) , X —>■ 0 

with L being slowly varying at 0. Hence it is possible to estimate a/A by 
applying Hill’s estimator [9] to the sample X^^,- ■ ■ , X~^. 


2.2 Nonparametric estimation 

Let ipn{u) = n~^ Y^=i denotes the empirical characteristic function (e.c.f.) 
obtained from the observations and its derivative. From (13), we are tempted 
to plug the e.c.f. of the observations to estimate the p.d.f. 9. Let (/i„)„>o 
and(/c„)n>o be two sequences of positive numbers such that 


lim hn = lim = 0 , 

n—^oo n—¥oo 


and consider the following sequence of estimators: 


9n{x) = max | ^ 


A ipn(u) 


\(pn.{u)\>Kr 


du, 0 


(14) 


Remark 2.2. We estimate l/ip{u) by l{|(,3„(u)|>K„}/‘i5n(u) with a suitable 
choice of a sequence (Kri)ra>i which converges to zero. The constant Kn is cho¬ 
sen such that \(pn{u) — '^(u)\ remains smaller than |:,5„(u)| and |:/5(u)| with high 
probability in order to avoid large errors when inverting (pn(u). In [16], the au¬ 
thors deal with the empirical characteristic function of i.i.d. random variables. 
In this case, the deviations of i/n(i^„(u) — ^p{u)) are bounded in probability, 
hence, they use as an estimator of l/(p(u). Here we 

truncate the interval of integration R by [—, where hn is a bandwidth 
parameter. This allows us to bound the estimation error —0 in sup norm. The 
deviation of ^Jn{if}n{u) — >f{u)) on depends on hn, see Theorem 4.1. 

The resulting is then taken slightly larger than 
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In order to evaluate the convergence rate of our estimator, we consider par¬ 
ticular smoothness classes for the density 9. Namely we define, for any positive 
constants K, L, m and s > 1/2, 

&{K, L, s, m) = is a density s.t. J \y\'^'^'^9{y) dy < K , (15) 

[ {1 + lup)'’ \J^9{u)\^ du < L^\ , 

JR I 

where TO denotes the Fourier transform of 0 

T9{u) = [ 9{y) e-'^“ dy , u G K . 


Hence L is an upper bound of the Sobolev semi-norm of 9. Note also that under 
(SN-l)-(SN-2), 9 belongs to Q{K,L,s,m) is equivalent to assuming that 


E 


1^01 


4+m 


<K and / {l + \u\‘^y\(pYoiu)\'^du<L'^ 


In the following, under assumptions (SN-l)-(SN-2), we use the notation Pg and 
Eg, where the subscript 9 added to the expectation and probability symbols 
indicates explicitly the dependence on the unknown density 9. The following 
result provides a bound of the risk P6((||6>—^n||oo > Mn) for well chosen sequences 
(hn), {tin) and (M„), which is uniform over the densities 0 G 0(/G, L, s, to). 

Theorem 2.1. Assume that the process X= given by (1) satisfies the 

assumptions /SN-1/-/SN-.2/ for some positive constants X and a. Let K,L,m 
be positive constants and s > Let CK,L,m,x/a the constant defined by 

CK,L,m,X/o. = exp ^ 

and C be a positive constant such that 0 < C < CK,L,m,x/a ■ Set 

hn = and Kn = C {l + , 

and define 9n by (If). Then, for n > 3, the density estimator 9n satisfies 


sup Eg 

9^0{K,L,s,m) 


9 - 9n 


CX) J 


< Mn-(2«-l)/(4^+2+4A/a) log(„)l/2 ^ 


(16) 


where M > 0 is a constant only depending on C,K,L,m,s and X/a. 

Proof. See Section 5.4. □ 


Remark 2.3. The constant C in Theorem 2.1 might be adaptively chosen. 
Indeed, the well known relationship (see [5] [Chapter 6] for example) between 
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the cumulant function of a filtered Poisson process and its intensity measure 
implies that the mean ^0 of Xq is given by 

fie = AE, [Fo] / = -Eg [Fq] < ^ 

Jo a a 

Since X is ergodic (see Section 4), the empirical mean /t„ of the sample 
Xi,--- ,Xn converges to fig almost surely and thus, replacing C in Theorem 
2.1 by Cn = exp(—/i„)/2, leads to the same rate of convergence since 

lim Cn = exp (- -Eg [Fq] ) /2 € {Q,CK,L,m,x/a) a.s. 

n->cx3 V •a / 

Remark 2.4. This theorem provides that the error in uniform norm converges 
at least at a polynomial rate 7^-(2 s-i)/(4s+2+8A/q!) log(n)^/^ that depends both 
on the quantity A/a and the smoothness coefficient s. For a given ratio A/a, 
the convergence rate becomes faster as s increases and tend to behave as 
when s —>■ 00 . On the other side, for a given smoothness parameter s, the 
rate of convergence decreases when the ratio A/a tends to infinity. This can 
be interpreted as the consequence of the pileup effect that occurs whenever the 
intensity A is large or the impulse response coefficient a is close to zero. 

Remark 2.5. Based on the previous theorem, one might wonder whether the 
rates of convergence are optimal. According to similar but not identical prob¬ 
lems ([16],[8]) in which authors estimate in a nonparametric fashion a Levy 
triplet (with finite activity) based on a low frequency sample of the associated 
process, the optimal rates of convergence are identical to ours. Our estimation 
procedure lies on stationary but dependent infinitely divisible random variables 
associated to an infinite Levy measure so that these results do not apply here. 
However we believe that the rates obtained in Theorem 2.1 are also optimal in 
this dependent context. The proof of this conjecture is left for future work. 


3 Experimental results 

The estimation procedure based on the estimator given by (14) can be made 
time-efficient and thus well suited to a very large dataset. In nuclear applica¬ 
tions, it is usual to deal with several million of observations while the intensity 
of the time-arrival point process can reach several thousand of occurrences per 
second. Typically, the shot-noise process in nuclear applications corresponding 
to the electric current is discretely observed for three minutes at a sampling rate 
of lOMhz and the mean number of arrivals between two observations lies be¬ 
tween 10 and 100. Such large values for the intensity and the number of sampled 
points motivate us to present a practical way to compute the estimator (14). 

Practical computation of the estimator 

In Section 2, we have defined the estimator of mark’s density by (14). Al¬ 
though it theoretically converges to the true density of shot-noise marks, the 


evaluation of the empirical characteristic function and its derivative based on 
observations Xi ,■■■ , Xn might be time-consuming when the sample size n is 
large. To circumvent this issue, we propose to compute the empirical character¬ 
istic function using the fast fourier transform of an appropriate histogram of the 
vector Xi, - ■ ■ More precisely, for a strictly positive fixed h, we consider 

the grid G = {hi : [minfc<„(Xfe)//iJ < I < |"maxfe<„(Xfe)//i]} and compute the 
normalized histogram H of the sample sequence (^i)i<;<n with respect to the 
grid G defined by 

1 " 

H{1) = - ^ 1[G(0;G(Z-Si)](^fe): Lmin(^fc)/^J < I < \max{Xk)/h'] - 1 . 


Denoting uin = [minfe<„(Xfc)//iJ and M„ = |"maxfe<„(Xfc)//i] — 1, remark 


that for every real u, we have: 

Ti -| I ^ j.yj. fi 

XI MG{iy,Gii+i)]{Xk) 


n Mn 




MXk 




k—1 l—rrin 


Replacing l[c{iy,G(i+i)]iXk) by l[G(i);G(i+i)] (^fc) foj. any real m, 

we get an approximation of the empirical characteristic function by defining 






Au) = Y. 


iu/i(Z+l/2) 


(18) 


l—r. 


For any real u, we have the following upper bounds 


\V>h,n{u) - 'y>n{u)\ < -|u| 


and 




Wh 


Xn{u) - ‘P'n{u)\ < ^ 1 + \u\h Y ^(0 + 1 / 2 ) 


l=m„ 


showing that the approximations are close to the true functions for small values 
of h and u. From these empirical characteristic functions, we construct an 
estimator of the marks’ characteristic function Lpy setting for any positive u: 

^^Xn{u) = 1 + ■ 

The advantage of using (fh,n is that ‘ph,n{u) ^'h ni'^) 1*® evaluated on a 

regular grid using the fast Fourier Transform algorithm. 

The last step in the numerical computation of the estimator (14) consists in 
evaluating the quantity 


'AY,h,niu)du = 




Using the Inverse fast Fourier Transform, we approximate the integral on a 
regular grid a; G by a Riemann sum. 
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Numerical results 


We now illustrate the finite sample behavior of our estimator on a simulated data 
set when the marks (in keV energy units) density follows a Gaussian mixture 
with 


p= [0.3 0.5 0.2] , /i=[4 12 22] , cr = [1 1 0.5] . 


Furthermore, in order to fit with nuclear science applications, where detectors 
have a time resolution of 10 Mhz, corresponding to a sampling time A = 10“^ 
seconds. The parameters of the experiment are set to a = 8.10® , A = 10®. 
Moreover, in nuclear spectrometry, the bandwidth hn is directly related to the 
known precision of the measuring instrument. For the following numerical ex¬ 
periment, we set it to 2.5 which is in range with the detector resolution as 
described in [14], Chapter 4. Figure 1 below shows a simulated sample path of 
such a shot-noise with its associated marked point process. 



Figure 1: Simulated Shot-Noise 


As shown in Figure 2 below, our estimator defined by (14) well retrieves 
the three modes of the Gaussian mixture as well as the corresponding variance 
from a sample of size 10®, which corresponds to a signal observed for one hun¬ 
dredth second. Current estimators used in nuclear spectrometry for similar data 
requires much longer measurements (up to 10 seconds). The reason is that these 
estimators do not consider observations where pile-up is suspected to occur, thus 
throwing away a large part of the available information. 
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Figure 2: Gaussian mixture case 


is provided in Table 1 


Moreover, an estimation of the risk Eg ||0 — 9n\\c 
for the shot-noise configuration described at the beginning of this section and 
for three different sample sizes. 


Sample size (n) 

Mean error (variance) 

1(F 

0.1015 (5,10.10-4) 

10^ 

0.0741 (9,51.10-^) 

To® 

0.0622 (9,38.10-'^) 


Table 1: Mean l°° error Eg 


11^* - ^nlloo 


for different sample sizes n with 9 given 


by the Gaussian mixture defined above. The displayed values are obtained by 
averaging ||0 — 9n\\oo over 100 Monte Garlo simulations. The variance of the 
error over the Monte Carlo runs is indicated between parentheses to assess the 
precision of the given estimates. 
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4 Error bounds for the empirical characteristic 
function and its derivatives 


To derive Theorem 2.1, since our constructed estimator involves the empirical 
characteristic function and its derivative, we rely on deviation bounds for 


Eg 

sup 



uG[—h~^,h~^] 



= 0,1 , h > 0 , (19) 


which are uniform over 9 G 0(K, L, s,m), where the smoothness class 
0(K, L, s,m} is defined by (15). Here, and ipn^ respectively denote the 
fc-th derivative of the characteristic function and its empirical counterpart as¬ 
sociated to the sample Xi ,..., These bounds are of independent interest 
and therefore are stated in this separate section. Upper bounds of the empirical 
characteristic function deviations have been derived in the case of i.i.d. samples: 
[8] [Theorem 2.2.] provides upper bounds of (19) for i.i.d. infinitely divisible ran¬ 
dom variables, based on general deviation bounds for the empirical process of 
i.i.d. samples found in [22]. Here we are concerned with a dependent sample 
Xi, • • • , Xn and we rely instead on [7]. We obtain upper bounds with the same 
rate of convergence as in the i.i.d. case but depending on the /3-mixing coeffi¬ 
cients, see Theorem 4.1. An additional difficulty in the non-parametric setting 
that we consider is to derive upper bounds that are uniform over smoothness 
classes for the density 9, and thus to carefully examine how the /3 coefficients 
depend on 9, see Theorem 4.2. 

Let us first recall the definition of /3-mixing coefficient (also called abso¬ 
lutely regular or completely regular coefficient) as introduced by Volkonskii and 
Rozanov [23]. For A,B two cr-algebras of ft, the coefficient I3{A,B) is dehned 
by 

/3(A, H) = i sup ^ |P(A, n B,) - F{A,)F{B,)\ 

the supremum being taken over all finite partitions and {Bj)j^j of H 

respectively included in A and B. When dealing with a stochastic process 
{Xt)t>o, the /3-mixing coefficient is defined for every positive s by: 


/3(s) = sup/3 {(t{Xu,u< t ), (ct {Xs+u,u > t)) 

t>o 


The process {Xt)t>o is said to be /3-mixing if limt_).oo /3(t) = 0 and exponentially 
/3-mixing if there exists a strictly positive number a such that /3(t) = 0(e““*) 
as t —>■ oo. 

We hrst state a result essentially following from [7] which specify how the 
/3 coefficients allows us to derive bounds on the estimation of the characteristic 
function and its derivatives. 


Theorem 4.1. Let k be a non-negative integer and Xi, • • • , X^ a sample of a 
stationary j3-mixing process. Suppose that there exists C > 1 and p G (0,1) such 
that (in < Up" for all n> 1. Let r > 1 and suppose that E [|Aii< oo. 
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Then there exists a constant A only depending on C, p and r such that for 
all h > 0 and n > 1, we have 


E 

sup 






max (e ,E [|Xi|2(fc+i)-] (l + ^log(l +/i-i)) 

<A-^- L. (20) 


Proof. The proof is deferred to Section 5.3 


□ 


It turns out that the stationary exponential shot-noise process X defined 
by (1) is exponentially /3-mixing if the first absolute moment of the marks is 
finite, see [15] [Theorem 4.3] for a slightly more general condition. However, in 
order to obtain a uniform bound of the risk of our estimator over a smoothness 
class, a more precise result is needed. In the sequel, we add a superscript 9 to the 
/3-mixing sequence to make explicit the dependence with respect to the mark’s 
density 0. The following theorem provides a geometric bound for the /3-mixing 
coefficients of the shot noise which is uniform over the class Q{K, L, s, m). 

Theorem 4.2. Let Xi, • • • , X„ be a sample of the stationary shot-noise process 
given by (1) satisfying (SN-3^-(SN-.8/ Let K,L,m > 0 and s > 1/2. Then 
there exist two constants C > 0 and p € (0,1) only depending on A, a, K, L, s 
and m such that, for all n > 1, 

sup /3® < Cp" < oo . (21) 

e&e{K,L,s,m) 

Proof. See Section 5.2. □ 


As a corollary of Theorems 4.1 and 4.2, we obtain error bounds for the 
empirical characteristic function when dealing with observations Xi, • • • , X„ of 
the stationary shot-noise process given by (1). 


Corollary 4.1. Let Xi, • • • , X„ be a sample of the stationary shot-noise process 
given by (1) satisfying /SN-1/-/SN-^/. Let K, L,m > 0 and s > 1/2 and let k 
be an integer such that Q < k < \ -\- n/2. Then there exists a constant B only 
depending on k, X, a, K, L, s and m such that for all and n > 1, 


sup Eg 

sup 


9ee{K,L,s,m) 




< B 


1 -b A/Ioi(TT7FT) 

ni/2 


This result can be compared to [8] [Theorem 2.2]. Note however that although 
our sample has infinitely divisible marginal distributions, it is not independent 
and the Levy measure is not integrable. 
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5 Proofs 

5.1 Preliminary results on the exponential shot noise 

We establish some geometric ergodicity results on the exponential shot noise 
that will be needed in other proofs. 

Definition 5.1 (Geometric drift condition). A Markov Kernel P satisfies a 
geometric drift condition (called DiV, /i, b)) if there exists a measurable function 
y : R —>■ [l,oo[ and constants {fJ,,b) € (0,1) x R+ such that 

PV <fiV + b. 

Definition 5.2 (Doeblin set). A set C is called a {rn,e)-Doeblin set if there 
exists a positive integer m , a positive number e and a probability measure v on 
M such that, for any x in C and A in ,B(M) 

P"^ix,A) > ev{A) . 

The following proposition is borrowed from [21] and relates explicitly the 
geometrical drift condition to the convergence in Wnorm (denoted by j • jy) to 
the stationary distribution. 

Proposition 5.1. Let P be a Markov kernel satisfying the drift condition 
D(y,p,,b). Assume moreover that for some d > 26(1 — /i) — l,mGN — {0} 
and e € (0,1), the level set {V < d} is an (m, e)-Doeblin set. Then P admits a 
unique invariant measure tt and P is V-geometrically ergodic, that is, for any 
0 < u < e/ {bm + fi^d — 1 + e)_|_ V 1 , n S N and x S R, 

||P"(x,-) - 7r||y < c{u)[k{V) + l/(x)]pL"/H(y) (22) 


where 

mb— 

• Om — minv i_^ 

• c(u) = u~^(l — It) + /r™ + 6m 

• p(u) = (1 - e + u(bm + p^d + e - 1)) V (l - 

• 7r(V) = ¥{Z £ V) where Z is a random variable with distribution tt 

In order to apply such a result to the sample [Xi, • • • , Xn) of the exponential 
shot-noise defined by (1), observe that it is the sample of a Markov chain which 
satisfies the autoregression equation = e~°‘Xi-\-Wij^i , where the sequence 

of innovations (Wi)igz is made up of i.i.d. random variables distributed as 

N>,{[0,1]) 

Wo= Y. . ( 23 ) 

k=l 


where 
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• -/Va([0, 1]) is a Poisson r.v. with mean A, 

• (^i)i>i i-i-d. r.v.’s with probability density function 0, 

• are i.i.d. and uniformly distributed on [0,1], 

• all these variables are independent. 

In the following, we denote by Qg the Markov kernel associated to the Markov 
chain under (SN-l)-(SN-2). 

Proposition 5.2 (Uniform Geometric drift condition). Let K,L,m > 0 and 
s > 1/2 and let 0 G Q{K, L, s, m). Then the Markov kernel Qg satisfies the drift 
condition D{V, 11 , b), where 


U : a: ^ 1 + |a;| , ^ = 6““ , 6=1 + - 6"“ . (24) 

Proof. We have for all 0 € &(K, L,s,m) and a: £ M, 

QgV{x) = Eg [l + |e““x + tUo|] 


< e-°‘V{x) + 1 - e-“ + = ^V{x) + b . 


□ 


Remark 5.1. A similar result holds for the functions U : a; —>■ 1 + |x/ where 
i G {I,--- , L4 + toJ}. 

Proposition 5.3 (Doeblin set). Let I > 1, K, L,m > 0 and s > 1/2 and define 
V as in (24). There exists e > 0 only depending on I, a, X, K, L,m > 0 and s 
such that, for all 0 G &{K, L, s,m), the Markov kernel Qg admits {V < 1} as 
an {1, e)-Doeblin set. 

Proof. Let 0 G Q(K, L, s,m). Denote by 0 the density of random variable 
with Ui and Yi two independent random variables respectively dis¬ 


tributed uniformly on [0,1] and with density 0. It is easy to show that, for all 
u G R, 



(25) 


The distribution of Wq is thus given by the infinite mixture 



where 6o is the Dirac point mass at 0 and 0*^ denote the fc-th self-convolution 
of 0. It follows that, for all Borel set A, 


Qeix,A)=e-^lAie-‘^x)+ / (1 - e-^)/e(e + e-“x) d^ 


A 
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In order to show that {V < Z} is a (e, l)-Doeblin-set for the kernels Qe, it is 
sufficient to exhibit a probability measure z/ such that, for all |a;| < 1 — 1 and all 
Borel set A 


fe{^ + e-“x)de > (1 - e-^)-ie z. {A) . 


Hence if for each 6 G 0(A", L, s, m) we find c(0) < d{9) such that 


inf inf inf \d(9) — ciO)] fg{^ + e “x) > 0 

ee0iK,L,s,m)c{e)<^<d{e)\x\<i-i 


the result follows by taking v with density [d{0) — c(0)]“^l[c(6i),d(e)] and e = 
(1 — e“'^)e'. By definition of fg above, it is now sufficient to show that there 
exist c{0) < d{0) and k{0) > 1 such that 

inf inf inf [d{0) - c{9)] ^ + e'^x) > 0 . 

ee0(K,L,s,m) c(e)<i<d(0) lxl<l-l 


Observe that for c < ^ < d and |x| < Z — 1 we have ^ + e““x G [c—e““(Z — 1), d + 
e““(Z — 1)]. So for any interval [c', d'] of length d' — c' > 2e““(Z — 1), we may 
set c = c' + e““(Z — 1) < d = d' — e““(Z — 1) so that c < ^ < d and |x| < Z — 1 
imply ^ + e““x G [c',d']. Hence the proof boils down to showing that for each 
0 G Q{K, L, s,m), there exist c'{0) < d'{0) with d'(d) — c'(d) > 2e““(Z — 1) and 
k{0) > 1 such that 


inf inf \d'(0) — c'(0) — 2e “ 

9ee{K,L,s,m) c'(0)<^<d'(0) 


(z-1)]>0. 


(27) 


By Lemma A.3, there exists a>0, A>0, d>l and eo > 0 and such that eg 
and A = b — a only depend on m, K, L and s (although a may depend on 0), 
and 

inf 0[x) > eo ■ 

a<.x< (a+A )/6 

Finally, Lemma A.4 and the previous bound yield (27), which concludes the 
proof. □ 


5.2 Proof of Theorem 4.2 

As explained in Section 5.1, {Xi)i>o is a stationary H-geometrically ergodic 
Markov chain with Markov kernel denoted by Qg. By [6], the /3-coefficient 
of the stationary Markov chain {Xi)i>o can be expressed for all n > 1 and 
0 G Q{K, L, s, m) as 


/3n = / IIQe(a^G) - TTflllTy 7re(dx) , 

Jr 

where irg is the invariant marginal distribution and | • \tv denotes the total 
variation norm, i.e. the H-norm with V = 1. Combining Propositions 5.2, 5.3 
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and 5.1, we can find constants C > 0 and p € (0,1) only depending on A, a, 
K, L, s and m such that 



(28) 


which concludes the proof. 

5.3 Proof of Theorem 4.1 

In [7], the authors establish a Donsker invariance principle for the process 

{Zn{f), / £ J^} where — P) is the normalized centered 

empirical process associated to a stationary sequence of /3-mixing random vari¬ 
ables {Xi, • • • , Xn) with marginal distribution P and is a class of functions 
satisfying an entropy condition. To be more precise, suppose that the sequence 
(X^) i>i is /3-mixing with ^ mixing rate function /3 is defined 

by /3(t) = /3pj if t > 1 and /3(t) = 0 otherwise while its cadlag inverse /3“^ is 
defined by: 


/3 ^{u) = inf{/3(t) < m} 


Further, for any complex-valued function /, denote by Q/ the quantile function 
of the r.v. |/(Xo)| and introduce the norm: 



The space £2,/? is defined as the class of functions / such that ||/||2,/3 < 00. In 
[7], the authors proved that (/I2./3, || • lb,/?) is a normed subspace of £2- A useful 
and trivial result from the definition of the norm £ 2^/3 (P) provides the following 
relation: 


I/I < I5K II/II2./3 < II5II2./3 . 


For any real r > I, another useful (less trivial) result in [7] states that under 
the condition 



we have £ 27 - C £ 2 ,/? with the additional inequality 



||/||2./3<||/||2. /l + r^/3„n-/(-i) 


n >0 


(29) 


where ||/||2r = E [|/(Xo)b’’] denote the usual 


-norm. 
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Now, we can state a result directly adapted from [7] [Theorem 3] that will 
serve our goal to prove Theorem 2.1. For the sake of self-consistency, we recall 
that, given a metric space of real-valued funtions (£1, || • ||) and two functions I, 
u , the bracket [/, u] represents the set of all functions / such that I < f < u. 
For any positive e, [Z,u] is called an e-bracket if ||1 — u|| < e. 

Theorem 5.1. Suppose that the sequence is exponentially /3-mixing 

and that there exists C > 1 and p G (0,1) such that /3„ < Cp^ for all n > 1. 
Let (T > 0 and let J- C be a class of functions such that for every f in J-, 
\\fh,p < <y. Define 

4>{a) = ^1-blog (iV[] (u, J', II • ||2,/3)) du , 

where || • || 2 ,/ 3 ) denotes the bracketing number, that is, the minimal 

number of u-brackets with respect to the norm || • || 2„0 that has to be used for 
covering T. Suppose that the two following assumptions hold. 

(DMRl) T has an envelope function F such that ||F'|l 2 r < oo for some 
r > 1. 


(DMR2) < oo . 

Then there exist a constant A > 0 only depending on C and p such that, for 
all integer n, we have 


E 


sup \Z, 


m 


Having this result at hand, we now remark that (19), for a fixed integer k, 
can be rewritten as 


E 

sup 


= n-i/^E 

sup |^„(/)| 







where 

= [-h-\h-^]}. (32) 

The proof of Theorem 4.1 based on an application of the previous theorem is as 
follows. 


Proof. We apply Theorem 5.1 for a fixed integer k, F = Fj/, F = Fk and 
r = (4 -b m)/4 where Fk : x ^ \x\^. 


Assumption (DMRl): Let A: be a fixed integer. On the one hand, the func¬ 
tion Fk is an envelope function of the class Ff and on the other hand, for any real 
r > 1, from (29), we have 


l|Ffc||2,/3 < E ^ /l+r^,0n n’'/!’’ b = < oo . 


(33) 
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Assumption (DMR2): For k a fixed integer, the class is Lipschitz in the 
index parameter: indeed, we have for every s,t in [—and every real x 


|(ia;)V""-(ia;)V‘"| <|s-t|| 


|fe + l 


(34) 


A direct application of [22] [Theorem 2.7.11] for the classes J-^ gives for any e > 0: 
A[] {2e\\F^+^\\2,p,Ft II • ||2,;3) < N (e, [-h-\h-^] , | • |) 

2/1-1 


< ! + ■ 


(35) 


where N and Ai[ j are respectively called the covering numbers and bracketing number 
(these numbers respectively represent the minimum number of balls and brackets of 
a given size necessary to cover a space with respect to a given norm). From (35), it 
follows that for any n > 0, we have 

= / A/l + log {u,F^,\\ ■ ||2,/3))du 

Jo 


< 


< 


A + log 1 + 


4||F’,+i||2,;3/r-i 


dtt 


(36) 




Jo \ ) 

= (7 + \\Fk+i\\y]ih~^^'^ < oo 


(37) 

because we supposed Fk+\ & C-ir and fin < Cp" which, from (29), implies that 
||fl)c+l||2,p < OO. 

Conclusion of the proof The application of Theorem 5.1 gives 

r 0(n'r,fc) 


Ee 


sup 




< A- 


, 1/2 


where A = A(1 + 1/Vl — )—i) since, from (33), we have ||Ffe||27-/o'r.,fc < 1. 
Set p From (36) and (33), we can write 


A + log 1 + 


<(>(o-) < 


For a = ar,k, we get after the change of variable v 


4| I I \2rC^^^h 


du , 


M\Fk + l\\2,rC^^pcr^f.k 


-du 


4>{’^r,k) < max (jj^feibr, ||F'fc + l|| 2 r)c^J 1^1 -b /l A/log(l -bu) 

By Lemma A.6, we get for a universal constant B > 0 that 

(l>{o-r,k) < Bmax (||Tfe||2r||Ffe + l||2r) ^1 -b \/log (1 + /i-i)^ . 

In the particular context of Corollary 4.1, we use the fact that ||Ffc|[ 27 - can be 
bounded by max(l, and ^ by a constant only depending on the param¬ 
eters K, L,s,m □ 
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5.4 Proof of Theorem 2.1 

Proof. We denote by 0° the function defined by : 


0O(x) ^ max ( 0, ^ / e-->y„(M)d^ 


1 




( 38 ) 


Since 9 G &{K, L, s, m) and s > 1/2, we have that P [d] is integrable and, under 
0 > 0, we have 

(39) 

(40) 


0(a:) = max (^0, — [ e (u)dz 

V 27r Jr 


We decompose the error in infinite norm as 

\\e-k\\^<\\9-e°j^ + \\9°-k\U 

From (38) and (39), we get 

\\0-On\\oo<- \ipYoiu)\du 


1 /■°° 

— \u ^u‘‘ipYo{u)\du 


< 


ImI ^®du 


1/2 


< 


L h' 


s-l/2 


\u‘‘lpyo{u)\^ du 


-(2s-l)/(4s+2+4A/a) 


1/2 


TT 2s — 1 (2s — l)7r ^ ^ 

where we used the Cauchy-Schwartz inequality and the assumption that 6 G 
&{K, L, s, m). 

We conclude with a bound of the term involving ||0° —0„||oo in (40). To this end, 
the following inequality will be useful. Using (12) and the mean-value theorem, 
we have 

P'xoi'^) 

~ O' „gR 


sup 

ttSR 


< - sup \ip'y^{u)\ < -E, [lUol] < . (42) 


V’Xo (u) 

By (12), we can bound the term ||0° — 0„||oo by 

Q 


\K-0n\\oo = 


Xtt 


PXoi^) Pn(u) 

- - ^\<P„{u)\>K„du 


< 


Xtt , 


2ah„^ 

< . sup 


PXoiu) Pniu) 
P'xoi-^) Pniu). 

^xoin) 

P'xoi'^) P'n{u) 


— h 


< 


A 


Xtt 


2ah^ 


\u\<h„ 


PXoiu) (finiu) 


du 


Att 


sup 

PXoi'^) P'niu) 

sup 

P'xoi'^) 


PXo[u) (fniu) 

\u\<hn^ 

PXq{u) 


\0n{u)\<K„ 


— 4n,l -|- An 
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Writing ^ as the term A^ i can be 

^Xoi^) 0niu) 


'£. 

V^Xq 

bounded as follows. 


sup 

\u\<h~ 

< < 


V’Xoiu) 0n[u) 


\(pn(u)\>Kn 


^ sup \^{u)\\ipn{u)-ipXo{u)\+K^^ 

\u\<hTi ^ 


sup \<P'n{'u) - ‘^'xoiu)\ 
\u\<hn ^ 


Thus, using (42), we get 


Eg [^n.l] < 


2/l„^Kn^Xl/(4+'")]Ee SUP|^|<^-l \(p'r,iu) - ^Xo(^)\ 


2ahj,^K ^Eg 


SUP|„|</,-i \<fn{u) - (pXo{u)\ 


ttA 


The two terms on the right hand side can be bounded using Corollary 4.1 with 
r = (4 + to)/ 4. It gives 


Eg 


sup \(p'n{'u) - ip'xo{u)\ 

\u\<h„^ 


< 


B if4/(4+m) ^ y^log(l + h-l)^ 


7.1/2 


and 


Eg 


sup \(pn{u) - ipXo{u)\ 

\u\<h^^ 


< 


b (^1 + ^log(l + /d^i)^ 


U/2 


In the following, for two positive quantities P and Q, possibly depending on 9 
and n we use the notation 


P ^ Q 4=^ for all n> 3, sup — < oo . 

eee(K,L,s,m) Q 


(43) 


(P is less than Q up to a multiplicative constant uniform over 6 £ 0(X, L, s, m)). 
We thus have that 


Ee [An,i 


< 


< 


1 + ^log(l + hn^) 
Kn n^l'^hn 
1 + ^log(l + hn^) 


„l/2;jV«+i 

< „-(2«-l)/(477+2+4A/a) ^ 


( 44 ) 


where we used the fact that <2C ^hn^°‘ for any integer n. 
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We now bound An^ 2 - From (42), remark that 



(45) 


From Lemma (A. 5), we have 


inf |(,5„(u)| > inf |v?jfo(M)|- sup \(pn{u) - (fXo{u)\ 
W\<hn \u\<hz'^ 



It follows that 


Pg inf \0n{u)\<Kn 
\ I’^l ^ 


<Pe[ sup \<Pn.{u)-ipXo{u)\>{CK,L,m,X/a-C){l + h^^) ) . 

\\u\<h-^ J 

Since 0 < C < Cx,L,m,\jaj applying Corollary 4.1 combined to the Markov’s 
inequality, and using (45), we get 



<^-(2s-l)/(4s+2+4A/a)iQg(^)l/2 ^ 


(46) 


Equations (41), (44) and (46) imply (16) and the proof is concluded. 


□ 


A Useful lemmas 

The following classical embedding will be useful. 

Lemma A.l (Sobolev embedding). Let K,L,m > 0 and s > 1/2. Let 0 S 
Q{K, L, s,m) defined in (15). Then, for any 7 S (0, (s — 1/2) A 1), there is a 
constant C > 0 depending on L, s and 7 such that, for every real numbers x, y. 


\eix)-0iy)\<C\x-yr , 


(47) 


where 



(48) 


The following result is used in the proof of Proposition 5.3. 
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Lemma A.2. Let K,L,m > 0 and s > 1/2. Let 7 G (0, (s — 1/2) A 1) and 
0 G 0(A, L, s, m). Then, there exists 0 < a < Tk such that 


inf e{x)>—{2K) 
a<a;<a+A lo 


l/(4+m) 


where Tk = ( 2 A)-i/ 0 ( 4 +m))^ ^ ^ C defined 

by (48). 

Proof. We first show that, for every T > 0, we have 

sup e{x) > ( 2 T)-i (1 - , 

|x|<T ^ 

Denote by D a random variable with p.d.f 0 belonging to the class Q{K, L, s, m). 
On the one side, we have 


P(|b^| <T)<2T sup 0{x) 

\x\<T 

and on the other side 

p (|r| < T) = 1 - p (|r| > T) > 1 - E [|y 1^+™] ^ 

where the first inequality is obtained via an application of the Markov inequality. 
Setting Tk = {2K f/'^^+'^\ we thus have 

sup 0{x) > {ATk)~^ . 

\x\<Tk 


Moreover, since 0 is continuous, we can without loss of generality suppose 
that there exists a positive number a in the interval (OjTr-] such that 

0{a) > (STx)-^ . 

From Lemma A.l, there exists a positive number A = (16TffC')“^/'^, indepen¬ 
dent of the choice of 0 such that 


inf 0{x) > (16Tr-) 
rcG[a,a+A] 


□ 


Lemma A.3. Let K,L,m,a > 0 and s > 1/2. Let 7 G (0, (s — 1/2) A 1) and 

0 G e{K,L,s,m). DefineTK = ( 2 A)-i/( 4 +™), A = 

with C defined by ( 48 ) and let 6 be a positive number satisfying 


1 < (5 < min 



Tk + A\ 

Tk J 


For any strictly positive v, define the function 0 by 0{v) = 
Then, there exists 0 < a < Tk such that 


inf 0(v) > 
a<'U<(a+A)/(5 


(2A)-i/(4+™)(5- 1) 
16a 


— 0{x)dx 

OiV Jv ^ ' 
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Proof. From Lemma A.2, we have 


inf e(x) > i-(2iF)-V(4+™) 
a<x<a+A 16 

for some a G (0,Tr:]. Let <5 G (l,e“ A ■ Since (a + A)/a is a decreasing 

function in a for a fixed A and 0 < a < Tk, we have that 


(a + A)/a > 


Tx + A 

Tk 


so that 5 < (a + A)/a . For any v G [a, (a + A)/5], we have 

pvd 


e{v) = 

> 

> 


A / 

av Jy 
vS — V 


CXV x^[v,v6] 

€k{S - 1 ) 


0{x)dx > 
inf 9{x) 


1 


av 


9{x)da 


which concludes the proof. 


□ 


The following elementary lemma generalizes the previous result for convolu¬ 
tions of lower bounded functions. 

Lemma A.4. Let 9, 9 two positive functions such that there exist positive num¬ 
bers a,b,c,d,e and i satisfying 

9{x) > el[Q_b](a;) and 9{x) > 

Then, for any 5 satisfying 0 < S < {b — a) A {d — c), we have 

{9*9) {x) > min(l,(5) ee l[a+c+s,b+d-S]{x) ■ (49) 

As a consequence, for any integer n in N*, we have 


^ (^) ^ ^min ^1, e '^[na-\-{b—a)/2,nb—{b—a)/2]{^') ■ 


(50) 


A lower bound of the decay of the absolute value of the shot-noise charac¬ 
teristic function is given by the following lemma. 

Lemma A.5. Assume that the process X is given by (1) under (SN-i)-(SN-^j 
with some positive constant a and A. Let K,L,m and s be positive constants. 
Then for all 9 G Q{K, L, s, m) and u G M, we have 


|‘/5Xo(w)| > CK,L,m,X/a (1 + kl) 


— A/c 


(51) 


where CK,L,m,\/a = exp(-A(L -K A:4/(4 +™))/c^). 
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Proof. From (9) and (11), we have for all u G K, 


«.(«)= exp (1/J^ 


e*’' - 1 


dr; I 9{x)dx I . 


If follows that 


bx.(»)| = exp(l/^(£ 


cos(r;) — 1 


dz;^ 9{x)dx^ 




First, we have for any real z and any function 9 G Q{K, L, s, m), 


|l-Re(v3vo(2))| = 


We thus get that 


[ (1 - e“^) 6»(a;)dx 

Js. 

< [ |l-e*“|6 »(a;)dx 

Js. 

<2 / |sin(a;z/2)| 0(a;)da; 

Js 

< f \xz\9{x)dx < 

JR 

Jo z 


Now, for |it| > 1, we have that 


< log |u| + z-l |Re (<^y„(z))| dz 

< log|u| +L , 


where we use the Cauchy-Schwartz inequality and 


|Re((/;yjr 


1/2 


< L . 


Inserting (53) and (54) in (52), we get the result. 


(52) 


(53) 


(54) 


□ 


Lemma A.6. There exists a constant B > 0 such that, for all u > 0, we have 
u i/log(1 + z;)-^ < B\/log (1 + u) 

J U ^ 
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Proof. For all m > 0, we have 


uf A/log(l +7;)^ = f y/\og{l+ uy)^ < Vu [ = 2\/u . 

Ju V ji y Ji y' 

As u —>■ 0, Vlog(l + u) is equivalent to ^/u. 

As M —?► oo, the Karamata’s Theorem (see [17] [Theorem 0.6]) applied to the 
function u —^ \/log(l + u)u~'^, which is regularly varying with index —2, gives 
that 

u / ^y\og{l + v)-^ ~ l/log (1 + u) , 
which concludes the proof. □ 
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